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Abstract. The combined effect of a dissipative fluid and quintessence 
energy can simultaneously drive an accelerated expansion phase at the 
present time and solve the coincidence problem of our current Universe. 
A solution compatible with the observed cosmic acceleration is succinctly 
presented. 



1. Introduction 

In recent years observation of Type la supernovae has lend strong support to 
an accelerated expansion of the Universe at present time Q. This unexpected 
feature could be explained by resorting to a small cosmological constant A, 
which on the one hand would provide enough negative pressure to account for 
this acceleration, and on the other hand would contribute an energy density 
of the same order of magnitude than the energy density of the matter content 
(baryonic plus dark) of today's Universe -say p m ~ 0.3 and p\ ~ 0.7. This 
seemly straightforward "solution" poses, however, a serious problem. Since p m 
redshifts as a~ 3 while A is a constant, why they turn out to be comparable 
today? i.e., why we happen to live in a very special (and rather short) phase of 
cosmic expansion? This is the so-called coincidence problem [Q]. 

To avoid it a new form of dark energy (quintessence energy, or Q-matter) 
has been introduced. This energy corresponds to a scalar field <p that slowly 
rolls down its potential with the key property of having a negative pressure. In 
some respects it mimics the scalar field that suppossedly drived inflation at the 
very early Universe -see e.g. [||. All these models -formulated for spatially flat 
Friedmann-Lemaitre-Robertson-Walker (FLRW) universes- overlook the fact 
that the matter component of the Universe (baryonic and non-baryonic) might 
not be very well approximated by a perfect fluid since, in general it should behave 
as a dissipative fluid and therefore it must have a non-equibrium pressure that 
might be non-negligible. 

In this short report we shall show that the combination of quintessence 
energy (an evolving scalar field with negative pressure) and a perfect matter 
fluid cannot simultaneously drive the current accelerated expansion and solve 
the coincidence problem. However, when the matter fluid is no longer assumed 
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perfect this difficulty disappears altogether for open and flat FLRW universes 
0. 

At this point it is fair to say, however, that attempts to solve the coincidence 
problem have been also made by using scalar-tensor theories of gravity rather 
than general relativity H. We will not deal with these here. 

2. Quintessence plus perfect fluid 

The stress-energy tensor of a perfect fluid (normal matter, i.e., baryonic and 
non-baryonic) plus a scalar field reads 

T a b = (Pm + P<f> + Pm + Pcl>)UaUb + {Pm + P^Qab {u a U a = -l). (1) 

The equation of state for the normal matter is p m = ( r y m — l)p m , with the 
baryotropic index lying in the range 1 < 7 m < 2, and 

P4> = \tf + v{4>) , P<t> = \tf - v{4>) ■ (2) 

A corresponding equation of state for the scalar field can be written as = 
(70 — l)p,£, or what it is the same 

where for non-negative potentials V(<j>) one has < 7^ < 2. However, the scalar 
field can be properly interpreted as quintessence if the restriction 7^ < 1 is met. 
The Einstein equations for any FLRW universe take the form 

n m + + n k = 1, (4) 



n = n(n-i)(3j -2)h , (5) 

where £1 = Q m + fi^, S7 m = p m /p c , Pfi/Pc with p c = 3H 2 the critical 

density, and = —k/(aH) 2 . As usual H = a/a denotes the Hubble factor, and 
k the spatial curvature index. 

Likewise the evolution equation for the scalar field <j> + SHcp + V' = can be 
recast as 

tl^ = [2 + (37 - 2) fl - 3 1<t> ] fl^H , (6) 
where 7 is the average baryotropic index defined by 

7SI = 7 m fJ m + 7^. (7) 

The combined measurements of the cosmic microwave background temper- 
ature fluctuations and the distribution of galaxies on large scales suggest a flat 
or nearly flat Universe ||. Hence the interesting solution of (||) at late times is 
SI = 1 (i.e., k = 0), and so we discard the solution f2 = as incompatible with 
observation. The solution S7 = 1 will be asymptotically stable for expanding 
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universes provided that condition 80,/dO, < holds in a neighborhood of f2 = 1, 
and this implies 7 < 2/3. Hence the matter stress violates the strong energy 
condition (SEC) p + 3p > and as a consequence the Universe accelerates its 
expansion, i.e., 3d /a = — (p + 3p)/2 > 0. 

Since the mixture of Q-matter and perfect dark matter fluid must violate 
the SEC, 70 must be low enough. Namely, because of 7 < 2/3, "y m > 1, and 
J(p < Ira-, equation (|7|) implies 7^ < 7. Then, introducing £1 = 1 in equation (||) 
we obtain 

^4, = 3(7 - ~f^H, (8) 

and therefore Cl^ > 0, i.e., will grow until the constraint (Q) is saturated, 
giving SIa = 1 in the asymptotic regime. Thus the matter fluid yields a van- 
ishing contribution to the energy density of the Universe at large times. This 
implies that a flat FLRW universe driven by a mixture of normal perfect fluid 
and quintessence matter cannot both drive an accelerated expansion and solve 
the coincidence problem. Therefore some other contribution must enter the 
stress-energy tensor of the cosmic fluid. A sensible choice is a negative pres- 
sure arising from the dissipative character of the matter component. It is worth 
mentioning that in deriving the above result neither j m nor 7^ were restricted 
to be constants. 



3. Quintessence plus dissipative fluid 

The only dissipative pressure that may enter the strees-energy tensor is a scalar 
pressure ir which has to be semi-negative definite for expanding fluids to comply 
with the second law of thermodynamics. In this light the stress-energy tensor 
keeps the same form as (|l]) but with p m replaced by p m + it. Now equation (|j) 
remains in place but 
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n = n(n - 1) 



37 + - 
p 



H 



(9) 



2 + 



7T 

3 7+~ 
P 



Q - 37^ > HQ^, 



(10) 



substitute equations (|5|) and (||) , respectively. The energy conservation equation 
of the normal matter is p m + 3 (7™ + ir/p m ) PmH = 0. Owing to the presence 
of the dissipative pressure it the constraint 7 < 2/3 does not longer have to be 
fulfilled for the solution O = 1 of equation (||) to be stable. Likewise, inspection 
of (|i~0|) shows that when Q = 1 one can have tl^ < just by choosing the ratio 
tt/p sufficiently negative. Thereby the constraint (||) allows a nonvanishing Q m 
at large times. By contrast tracker fields based models (valid only when 0,^ = 0) 
predict that Q m — > asymptotically. 

A fixed point solution of equation (g) is O = 1. Note that equations (||) and 
( |l0| ) have fixed point solutions Q m = Q m o and = fJ^O; respectively, when the 
partial baryotropic indices and the dissipative pressure are related by 
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tt 2H 

Then, the smaller 7^, the larger the dissipative effects. Let us investigate the 
requirements imposed by the stability of these solutions. From (|9|) we see that 
7 + 7T//? < 2/3 must be fulfilled if the solution 0, = 1 is to be asymptotically 
stable. This condition, together with (|TT]) , leads to the additional constraint 
on the viscosity pressure tt < (2/3 — j m ) p m , which must be negative. Also by 
virtue of (||) and the first equality in (|Tl| ) we obtain from last relationship that 
7^<2/3. 

In the special case of a spatially flat universe ($7 = 1), the stability of the 
solutions J7 m o and Q^q may be studied directly from (|To|). Setting = Q^ + lj, 
with I u |<C f^o, and using (^) it follows that 
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3^ m ( 7m -7 + — ) H (n^o + to) . (12) 



Accordingly the solution Q = 1, $7^ = Q^q is stable for the class of models that 
satisfies ip = 7 m — 7^ + (7r/p m ) < and -0 —> for t — > 00. Note that this 
coincides with the attractor condition (|i~l|). 

To study the stability of the solutions f2 m o and S7^o when k 7^ we introduce 
the parameter e = fi m /fi^,. As it turns out its evolution is governed by 
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and perturbating this expression about the solution e = eo ~ 0(1), (i.e., using 
the ansatz e = €q + 5 with \5\ ^ 1) we obtain with the help of ( |ll|) 

* = -^-(|-7«)n*^(eo + «) (14) 

near the attractor. For 0,^ > (negatively spatially curved universes) the ratio 
(17 m /ri^)o is a stable solution. For < one has to go beyond the linear 
perturbative regime and/or restrict the class of models as in the spatially flat 
case to determine the stability of the solution. 

We would like to stress that by large times we mean times after the cos- 
mological perturbations evolved into the nonlinear regime. Thus the structure 
formation scenario will not be spoiled by the quintessence field. 

Recently there have been claims that CDM should not be a perfect fluid 
because it ought to self-interact (with a mean free-path in the range 1 kpc < 
I < 1 Mpc) if one wishes to explain the structure of the halos of galaxies 0]. 
It is not unreasonable to think that this same interaction lies at the root of the 
dissipative pressure tt at cosmological scales. Bearing in mind that / = 1/na, 
with n the number density of CDM particles and a the interaction cross section, 
a simple estimation reveals that at such scales I is lower than the Hubble distance 
H^ 1 and accordingly the fluid approximation we are using is valid. 

There exist a handful of solutions with the desired asymptotic properties 
-see [|| for details. By way of example we just mention a ~ i~ 2 / A with 
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A = ~ j- (37 + u) + [(3 7 - vf + 36 7m v 2 n m ] 1/2 j , (15) 

where v denotes the number of interactions between CDM particles in a Hubble's 
time and v the speed of the dissipative signals. 

We may conclude by stressing that both acceleration and coincidence can 
be satisfactory explained by a combination of quintessence and dissipative dark 
matter. For these models attractor solutions exist with very interesting proper- 
ties: an accelerated expansion, spatially flatness and a fixed ratio of quintessence 
to dark matter energy density. In consequence, the quintessence scenario be- 
comes more robust when the dissipative effect of the nonequilibrium pressure 
arising in the CDM fluid is allowed into the picture. 

In Q we presented specific models with an ample region in the space of out- 
of-equilibrium thermodynamic parameters satisfying observational constraints 
in the asymptotic attractor regime which our Universe may well be approaching. 
In a future research we shall aim to improve these constraints by analysis of the 
luminosity distance-redshift relation for type la supernovae and simulations of 
structure formation that include dissipative effects. 
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